Isotropic-nematic transition in hard-rod fluids: relation between continuous and restricted-orientation models

We explore models of hard-rod fluids with a finite number of allowed orientations, and construct their bulk phase diagrams within Onsager's second virial theory. For a one-component fluid, we show that the discretization of the orientations leads to the existence of an artificial (almost) perfectly aligned nematic phase, which coexists with the (physical) nematic phase if the number of orientations is sufficiently large, or with the isotropic phase if the number of orientations is small. Its appearance correlates with the accuracy of sampling the nematic orientation distribution within its typical opening angle. For a binary mixture this artificial phase also exists, and a much larger number of orientations is required to shift it to such high densities that it does not interfere with the physical part of the phase diagram.Comment: 4 pages, 2 figures, submitted to PR

Free Minimization of the Fundamental Measure Theory Functional: Freezing of Parallel Hard Squares and Cubes

Due to remarkable advances in colloid synthesis techniques, systems of squares and cubes, once an academic abstraction for theorists and simulators, are nowadays an experimental reality. By means of a free minimization of the free-energy functional, we apply Fundamental Measure Theory to analyze the phase behavior of parallel hard squares and hard cubes. We compare our results with those obtained by the traditional approach based on the Gaussian parameterization, finding small deviations and good overall agreement between the two methods. For hard squares our predictions feature at intermediate packing fraction a smectic phase, which is however expected to be unstable due to thermal fluctuations. This implies that for hard squares the theory predicts either a vacancy-rich second-order transition or a vacancy-poor weakly first-order phase transition at higher density. In accordance with previous studies, a second-order transition with a high vacancy concentration is predicted for hard cubes

Charge regulation and ionic screening of patchy surfaces

The properties of surfaces with charge-regulated patches are studied using non-linear Poisson-Boltzmann theory. Using a mode expansion to solve the non-linear problem efficiently, we reveal the charging behaviour of Debye-length sized patches. We find that patches charge up to higher charge densities if their size is relatively small and if the patches are well separated. The numerical results are used to construct a basic analytical model which predicts the average surface charge density on surfaces with patchy chargeable groups.Comment: 9 figure

Phase behavior and interfacial properties of nonadditive mixtures of Onsager rods

Within a second virial theory, we study bulk phase diagrams as well as the free planar isotropic-nematic interface of binary mixtures of nonadditive thin and thick hard rods. For species of the same type the excluded volume is determined only by the dimensions of the particles, whereas for dissimilar ones it is taken to be larger or smaller than that, giving rise to a nonadditivity that can be positive or negative. We argue that such a nonadditivity can result from modelling of soft interactions as effective hard-core interactions. The nonadditivity enhances or reduces the fractionation at isotropic-nematic ($IN$) coexistence and may induce or suppress a demixing of the high-density nematic phase into two nematic phases of different composition ($N_1$ and $N_2$), depending on whether the nonadditivity is positive or negative. The interfacial tension between co-existing isotropic and nematic phases show an increase with increasing fractionation at the $IN$ interface, and complete wetting of the $IN_2$ interface by the $N_1$ phase upon approach of the triple point coexistence. In all explored cases bulk and interfacial properties of the nonadditive mixtures exhibit a striking and quite unexpected similarity with the properties of additive mixtures of different diameter ratio.Comment: 12 pages, revised version, submitted to JC

Density Functional Theory for Chiral Nematic Liquid Crystals

Even though chiral nematic phases were the first liquid crystals experimentally observed more than a century ago, the origin of the thermodynamic stability of cholesteric states is still unclear. In this Letter we address the problem by means of a novel density functional theory for the equilibrium pitch of chiral particles. When applied to right-handed hard helices, our theory predicts an entropy-driven cholesteric phase, which can be either right- or left-handed, depending not only on the particle shape but also on the thermodynamic state. We explain the origin of the chiral ordering as an interplay between local nematic alignment and excluded-volume differences between left- and right-handed particle pairs

Dynamical Heterogeneities and Cooperative Motion in Smectic Liquid Crystals

Using simulations of hard rods in smectic-A states, we find non-gaussian diffusion and heterogeneous dynamics due to the equilibrium periodic smectic density profiles, which give rise to permanent barriers for layer-to-layer diffusion. This relaxation behavior is surprisingly similar to that of non-equilibrium supercooled liquids, although there the particles are trapped in transient (instead of permanent) cages. Interestingly, we also find stringlike clusters of up to 10 inter-layer rods exhibiting dynamic cooperativity in this equilibrium state.Comment: 10 pages, 4 figure

The sediment of mixtures of charged colloids: segregation and inhomogeneous electric fields

We theoretically study sedimentation-diffusion equilibrium of dilute binary, ternary, and polydisperse mixtures of colloidal particles with different buoyant masses and/or charges. We focus on the low-salt regime, where the entropy of the screening ions drives spontaneous charge separation and the formation of an inhomogeneous macroscopic electric field. The resulting electric force lifts the colloids against gravity, yielding highly nonbarometric and even nonmonotonic colloidal density profiles. The most profound effect is the phenomenon of segregation into layers of colloids with equal mass-per-charge, including the possibility that heavy colloidal species float onto lighter ones

Depletion-induced biaxial nematic states of boardlike particles

With the aim of investigating the stability conditions of biaxial nematic liquid crystals, we study the effect of adding a non-adsorbing ideal depletant on the phase behavior of colloidal hard boardlike particles. We take into account the presence of the depletant by introducing an effective depletion attraction between a pair of boardlike particles. At fixed depletant fugacity, the stable liquid crystal phase is determined through a mean-field theory with restricted orientations. Interestingly, we predict that for slightly elongated boardlike particles a critical depletant density exists, where the system undergoes a direct transition from an isotropic liquid to a biaxial nematic phase. As a consequence, by tuning the depletant density, an easy experimental control parameter, one can stabilize states of high biaxial nematic order even when these states are unstable for pure systems of boardlike particles

Coupled water, charge and salt transport in heterogeneous nano-fluidic systems

We theoretically study the electrokinetic transport properties of nano-fluidic devices under the influence of a pressure, voltage or salinity gradient. On a microscopic level the behaviour of the device is quantified by the Onsager matrix ${\bf L}$, a generalised conductivity matrix relating the local driving forces and the induced volume, charge and salt flux. Extending ${\bf L}$ from a local to a global linear-response relation is trivial for homogeneous electrokinetic systems, but in this manuscript we derive a generalised conductivity matrix ${\bf G}$ from ${\bf L}$ that applies also to heterogeneous electrokinetic systems. This extension is especially important in the case of an imposed salinity gradient, which gives necessarily rise to heterogeneous devices. Within this formalism we can also incorporate a heterogeneous surface charge due to, for instance, a charge regulating boundary condition, which we show to have a significant impact on the resulting fluxes. The predictions of the Poisson-Nernst-Planck-Stokes theory show good agreement with exact solutions of the governing equations determined using the Finite Element Method under a wide variety of parameters. Having established the validity of the theory, it provides an accessible method to analyse electrokinetic systems in general without the need of extensive numerical methods. As an example, we analyse a Reverse Electrodialysis "blue energy" system, and analyse how the many parameters that characterise such a system affect the generated electrical power and efficiency

Sedimentation of binary mixtures of like- and oppositely charged colloids: the primitive model or effective pair potentials?

We study sedimentation equilibrium of low-salt suspensions of binary mixtures of charged colloids, both by Monte Carlo simulations of an effective colloids-only system and by Poisson-Boltzmann theory of a colloid-ion mixture. We show that the theoretically predicted lifting and layering effect, which involves the entropy of the screening ions and a spontaneous macroscopic electric field [J. Zwanikken and R. van Roij, Europhys. Lett. {\bf 71}, 480 (2005)], can also be understood on the basis of an effective colloid-only system with pairwise screened-Coulomb interactions. We consider, by theory and by simulation, both repelling like-charged colloids and attracting oppositely charged colloids, and we find a re-entrant lifting and layering phenomenon when the charge ratio of the colloids varies from large positive through zero to large negative values